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Abstract 

We develop the general formalism for the study of neutrino propagation in pres- 
ence of stochastic media. This formalism allows the systematic derivation of evolution 
equations for averaged quantities as survival probabilities and higher order distribu- 
tion moments. The formalism applies equally to any finite dimensional Schroedinger 
equation in presence of a stochastic external field. New integro-differential equations 
valid for finite correlated processes are obtained for the first time. For the particular 
case of exponentially correlated processes a second order ordinary equation is obtained. 
As a consequence, the Redfield equation valid for Gaussian delta-correlated noise is 
rederived in a simple way. The formalism, together with the quantum correlation 
theorem is applied to the computation of higher moments and correlation functions 
of practical interest in forthcoming high precision neutrino experiments. It is shown 
that equal and not equal time correlators follow similar differential equations. 



PACS numbers: 14.60.Pq; 13.10.+q; 13.15.+g; 14.60.Gh; 96.60.Kx; 02.50.Ey; 05.40.+j; 
95.30.Cq; 98.80.Cq. 
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1 Introduction. 



Neutrino oscillations in the presence of matter and magnetic fields have been an area 
of intense interest for a long time. In many astrophysical situations the matter density 
and/or magnetic field may fluctuate around a mean value. In some scenarios it is expected 
that these fluctuations are very strong, for example when dealing with the magnetic fields 
that the neutrino encounters while propagating in the convective regions of the Sun. 

The case of neutrino spin precession in a noisy magnetic field with a <5-correlated 
Gaussian distribution was considered initially by fl]] for neutrinos in vacuum. In a 
differential equation (a Redfield equation) for the averaged matrix density was derived for 
the case in which the random noise was taken to be again a delta-correlated Gaussian 
distribution. Numerous works where approximate or numerical approaches are used to 
study the implications of matter random perturbations upon the SNP problem, Supernova 
dynamics or others are avalaible (among them: [||, ||, ||, [7|, ||, [R]] see introductions in 

[7| for a more extensive review). 

In this work we will develop the general formalism for the study of neutrino propaga- 
tion in presence of stochastic media. This formalism will allow the systematic derivation 
of evolution equations for averaged quantities as survival probabilities and higher order 
moments of them. New integro-differential equations valid for finite correlation processes 
are obtained for the first time. For exponentially correlated processes a second order or- 
dinary equation is obtained as a consequence. The Redfield equation valid for Gaussian 
delta-correlated noise is rederived in a simple way as a particular case: it is obtained as 
the zero-order term of a asymptotic expansion in the inverse of the correlation length. In 
the context of neutrino oscillations, it is shown that the presence of matter noise induces 
the appearance of a effective complex part in the matter density. The presence of mag- 
netic field noise affects however in a qualitatively different, simpler, way. The range of 
validity of the 5-correlated model is checked in a realistic example corresponding neutrino 
oscillations in random solar matter. Finally the formalism, together with the quantum cor- 
relation theorem will be applied to higher moments and correlation functions of practical 
interest specially in forthcoming high precision neutrino experiments. 

The formalism may be applied to any quantum system governed by similar Schroedinger 
equations: equations where a stochastic function appears multiplicatively in some part of 
the Hamiltonian. For simplicity of notation, in this work we will always deal with Hamil- 
tonians which are finite dimensional operators in some Hilbert space, the equations will be 
obviously of the same type for general systems described by infinite dimensional Hamilto- 
nians. 

2 The General Case. 

Let's consider a general system whose evolution is described by the linear random Schroedinger 
equation 

id t X = p{t)L{t)X- X(0)=X (1) 
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where: X is a vector (i.e. a multiflavour wave function) or a matrix (i.e. a density matrix) 
of arbitrary dimension. Eq.(|l|) is defined in a interaction representation, any additive non- 
stochastic term has been solved for and absorbed in the definition of X(t) and L(t). L(t) 
is a general linear operator. It can be an ordinary matrix if X represents a wave function. 
If X is a density matrix, L(t) is a commutator: 

L A (t)X = [A(t),X] (2) 

Any other linear operator is admissible: obviously it is always possible to define an enlarged 
vector space where the action of L is represented by a matrix. 

We will suppose that p(t) in Eq.(jl]) is a scalar Gaussian process completely determined 
by its first two moments, which, without loss of generality can be taken as: 

(p(t)>=0; (p(t)p(t')) = f(t,t') (3) 

In other terms, p(t) is characterized completely by the measure 

1 r°° 

[dp]exp--/ p{t)r l {t,t')p(t')dtdt'. (4) 

An important particular case is when the process is <5-correlated, the correlation func- 
tion is of the form: 

f(t,t') = n 2 6(t-t'). (5) 

A convenient and practical way of parameterizing a correlation function with a finite 
correlation length is to use an exponential function: 

/(MO = ^ 2 eexp(-e | t-t' |) (6) 

In this case the correlation length is defined as r = 1/e. The expression (j|) is reobtained 
letting r — > or e — > oo in Eq.(j||). 

X(t), the solution to Eq.(@) is a stochastic function. The objective of this work is to 
obtain in a systematic way equations for its ensemble average and higher moments. Let 
us remark that for us Eq.([l|) is purely phenomenological, we suppose that is the result of 
a more complete microscopic analysis which can account for the randomness of p(t) (see 
for example 0). 

To obtain a differential equation for the average of X we make use of the following well 
known property: For any Gaussian process p(t) characterized by a 5-correlation function 
(pi) and any functional F[p] we have the following relation ([11]): 



For a Gaussian process with arbitrary correlation function (Eq.(|4|)) we have the general 
relation: 

(F[p]p(t))= Jdr(pmr))(f^). (8) 
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Let's consider now the evolution operator U for a particular realization of Eq.(Q); by 
definition: 

X(t) = U(t,t )X . (9) 

The operator U(t, to) is a functional of p and it has the following formal expression in 
terms of a time ordered exponential: 

U(t,t ) =Texp-i [ p(t)L(t)cIt. (10) 
Jto 

The functional derivative of U with respect to p(t) can be computed by direct methods 
differentiating term by term the series expansion for Eq.(|i~0|). The result is: 

S J^l = -iL(T)U(T,t ), t <T<t. (11) 

In order to obtain a differential equation for (X) we observe that 

i(d t X) = id t (X) = L(t)(p(t)X). (12) 



Combining together Eq. (p|,p|,pTD we obtain easily that the evolution equation for the en- 
semble average is in the general case a integro-differential equation given by: 

id t (X(t)) = -i tdt' f(t,t')L(t)L(t>)(X(t>)), (X(0))=X o (13) 



o 

This equation, which is exact and of very general validity, is the equation we were looking 
for and one of the main result of the present work. Note that it has been obtained 
previously integro-differential equations, valid for particular cases or in particular limits, 
using heuristic ad-hoc arguments (for example in |l(|). The derivation of Eq.(|i~3|) done 
here is the rigorous justification for such approaches. 



3 The particular (5-correlated case. 

For the particular case where the correlation function is of exponential type, a second 
order ordinary differential equation can be obtained as we will see below. In the other 
hand, for the simpler 5-correlated case the evolution equation ( |i~3| ) becomes the purely 
differential equation: 

id t (x(t)) = -m 2 L 2 (t)(x(t)). (14) 

Taking L as a commutator, this last equation coincides with the Redfield equation derived 
by Note that the effective " Hamiltonian" appearing in the second part of Eq.(|l4|) is 
not anymore hermitic ( an example of a fluctuation-dissipation effect). 

In practical cases of interest for the neutrino oscillation problem, the original equation 
for the density matrix is of the slightly simpler form: 

id t X = [H Q {t) + P {t)g{t)H u X] ; X{0) = X (15) 
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where we have written explicitly the commutator, pit) is a stochastic function as before 
and g{t) an arbitrary scalar function. Hq,Hi are Hamiltonian matrices, the former con- 
tains the average part of pit): H = H' (t) + po(t)Hi, with p = (p). The latter is supposed 
to be time independent. 

For the problem described by Eq.([l5|) the corresponding Redfield equation for the 
averaged density matrix is of the form 

id t (X) = [Ho(t),(X)]-in 2 g 2 (t)[Hi,[H h (X)]]. (16) 
= H (X)-(X)H+ + 2in 2 g 2 it)H 1 (X)H 1 (17) 

where in the last line the following effective Hamiltonians were defined: 

H± = H ±ig 2 it)H 2 . 

The solution of what is called the " coherent" part of Eq. ( |l7|) is accomplished by defining 
the average evolution operator: 

(U±) = TexpJdrH^ir); (IT) = (f/+) t . (18) 

The coherent part of the density matrix is then: 

(X)coh = {U-)X (U + y. (19) 
Defining a new "coherent" interaction representation by the relations: 

H P = {U-y^U-), H Q = (U-)H 1 {U-)-\ (X) I = (U-)(X) coh {U-)-\i20) 
the resolution of the original equation is equivalent to the resolution of the following one: 

id t (X)j = 2ig 2 it)n 2 H P (X) I H Q . (21) 

There are some important particular practical cases where Eq.(|l4D can be solved or 
simplified considerably taking into account the algebraic properties of a specific L (in what 
follows k(t) is always a scalar function): 

A. Let us suppose that L is such that L 2 it) = k(t)L(t). This case appears in the 
computation of the average wave function with matter density noise. Eq.(|i~4|) reduces to: 

id t (x) = -m 2 k(t)Lit)(x). 

The averaged equation is similar to the original one, the non-random part of the density 
is "renormalized" acquiring an imaginary term: 

P ~ * Po — iQ 2 k. 

This is the density which will appear in the coherent effective Hamiltonians Hq. This 
case was discussed already in B. Enquist and Semikoz (H) saw numerically that the net 
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effect in the averaged survival probability of the same system is a precession reduction, 
similar to the effect of a larger constant matter density. It can be argued (Eqs. (|l9| ) 
through (21)) that the complex renormalization defined before appears also, but this time 



approximately, if we compute average probabilities. Under this perspective the reduction 
in neutrino precession can be explained as a smoothening of the MSW resonance induced 
by the imaginary term. It is expected from here that (^-correlated matter noise have 
importance only if applied over a MSW resonance region. 

B. The case where L 2 (t) = k(t)I with / the identity matrix appears in the computation 
of the averaged neutrino wave function under noisy magnetic spin-flavour precession. The 
resulting equation can be integrated trivially (to be compared with the previous case): 

(X(t)) = exp-ft 2 f drk{r) X . 
Jo 

The average wave oscillation is damped by a factor equivalent to the one first calculated 
by Nicolaidis ( (!]] ) . This damping manifest itself also when computing the average density 
matrix from Eq.(|2~l|). From these difference of behavior with respect case (A) it is expected 
that magnetic field noise can affect even if applied far from any resonance region. 

C. The case where L A (t) = —k(t)L 2 {t) appears in the computation of full average 
density matrix with both, matter density or magnetic noise. We can obtain in this case 
the "conservation law": 

(l-tt 2 k{t)) L 2 (t)dt(X{t)} =0. 

L 2 (t) is not invertible because the operator L(t) has a zero eigenvalue. The previous ex- 
pression has proved to be of practical importance in some concrete numerical applications 
(0,11,0]). 

4 An Asymptotic Expansion for Exponentially correlated 
systems. 

We will see now how the Eq,(|i~4]) can be obtained as a limiting particular case when 
the correlation length tends to zero. For this purpose we use an exponential correlation 
function as Eq.(||), the integro-differential evolution equation becomes in this case: 

id t (X(t)) = -iQ 2 eexp(-et) f l dt'exp(eT)L(t)L(t')(X(t')). (22) 

Jo 



Let us compute the asymptotic expansion of the second term of Eq.(|22|) valid for e 
large, the following expansion is valid for any function g{t) 

h(e) = e exp(-et) f dr exp(er)g(r) ~ g(t) - ^ + ^ + . . . (23) 
Jo e e z 

Inserting this expression in Eq.(|22|), we obtain the following expansion in powers of e: 

id t (X) = -m 2 L 2 (t)(X) + i^L(t)d t (L(t)(X}) + o(±] (24) 
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To leading order in 1/e, we recover the expression corresponding to the 5-correlated case. 
At next-to-leading order we get finite-correlation correction terms: 

id t (X) = -m 2 L 2 {t)(X) +i^-L(t)L'(t)(X) + i^-L 2 (t)d t (X) (25) 
or equivalently 

- ^ 2 (*)) d t (X) = (-tfL 2 (t) + ¥-L(t)L'(t)J (X) (26) 

We finally get the following differential equation valid to order 1/e, making the sup- 
position that the operator which multiplies the first term is invertible, 

d t (X) = ^-n 2 L 2 (t) + ^L(t)L'(t) + ^-L\t)j(X) (27) 

This equation can be used for finite, but relatively large, correlation lengths. We see 
that, up to this degree of approximation, not only the ratio level of noise to correlation 
length (O 2 /e) is important. We have different regimes according to the value of f2 2 . For 
low noise amplitude (Q 2 << 1) the first term will be more important. For strong noise 
(f2 2 >> 1) the second one, proportional in this case to L 4 , will dominate. 

We have pointed out previously that in the practical cases of interest for neutrino 
oscillation problems the relation L 4 (t) = k(t)L 2 (t) holds, with k(t) positive. If in addition 
the term L'(t) can be neglected, as indeed happens in some occasions, Eq.(p7|) becomes: 

d t (X) = -fi 2 ^l - ^fc(t)) L 2 (t) (X) (28) 

In this case and within this level of approximation the net effect of the presence of a finite 
correlation length is visible: it amounts to a rescaling and consequent reduction of the 
noise parameter Q 2 . 

5 An Exact Differential Equation for Exponentially Corre- 
lated Systems. 

In contrast with the approximate approach used in the previous section, we can derive 
in fact a simple, ordinary second order differential equation for the case of exponential 
correlation. Let's suppose that the original equation is of the same decomposable type as 
the one appearing in Eq.(|l^) but let's include other cases using the general notation: 

id t X = (LoW+pflgW^X; X(0) = X (29) 

with Lq,L\ general linear operators as before, the latter time-independent. In this case 
Eq.([D|) is of the form 

id t (X) = L (t)(X)-in 2 eL 2 e^ t g(t) f ' dt 'g(t')e et '{X(t')). (30) 



o 



7 



Differentiating the equation once and after some simple algebra we obtain the following 
ordinary second order differential equation: 



(d t -\(t))[id t -L (t)]{X) = -iQ 2 eg 2 (t)L 2 {X), 

i(d t (X)) = L(0)(X) , (X) =X ; (31) 

with \(t) = — e + g'(t)/g(t). Let's remark that this equation is exact, unless the approxi- 
mations we have seen in the previous section. 

As illustration and for the sake of comparation, we have solved Eqs.(j3l|) and (|l4| ) for 
the case of two-flavour ultrarelativistic neutrino oscillations in presence of solar matter 
with a random perturbation. The function g(t) is in this case proportional to the solar 
electron density profile and can be approximated rather well by an exponential function 
(see notation and details in 0] for example). The results of the numerical calculations 
are shown in Figs.(^) (A) (Eq.(|l4|)) and (B) (Eq.(|3l|)). In plot (A) we observe the typical 
structure of the solar MSW resonance as a function of the neutrino energy and how this 
resonance is modified by the presence of the 5 correlated chaotic field. The broad range 
of validity of Eq.(14), considered as an approximation to the more realistic Eq.(31) is 



apparent in plot (B). We observe there how the survival probability depends on e. The 
neutrino energy is fixed in this case to be at the center of the resonance region and 
corresponds approximately to a neutrino oscillation length of 3.1 10 3 Km, 4 10 -3 Rq or 
8 10 -3 times the distance over the which the noise is acting. For any fixed value of the 
quantity Q 2 , the limit of small e (large L/d in the figure) coincide with the noise- less 
case. After a relatively narrow transition region, for larger, but in fact not so large, e ( 



L ~ d) the probability tends to a constant; the one which is obtained from Eq.(14). The 
conclusion is that, at least for the example in consideration, the 5-correlated noise is a 
good model even for "not so small" correlation lengths. [] 



6 Higher Order Moments and the Quantum Regression The- 
orem. 

Second order moments, expressions of the type (XiXj), or in general, moments of any 
order, can be computed also using equations similar to Eq.(li~3|), The straightforward 
procedure is to define products Xij.^ = XiXj.-X^ and write differential equations for 
them using the constitutive equations for each of the X{. The resulting equations are of 
the same type as Eq. ([[]). The similarity is obvious when one adopts a tensorial notation 
and defines products of the form X^ n ' = X <S> ... <S> X^ . The derivation of equations is 
specially simple within this notation. 

In the computation of averages of quantities of physical interest, i.e. expected signal 
rates, appear also correlators of quantities at different times. In the most simple case 

1 These results should be compared with the "cell" calculations presented in Jl5| . In such approach 
the averaged probability is computed by direct "Monte-Carlo" methods. In order to simulate a finite 
correlation length a step function correlation is used: f(t, t') is negligible whenever | t — t' | is greater than 
some characteristic scale, otherwise is a constant. 
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expressions of the type {Xi(t + r)Xj(t)). We will show shortly that the knowledge of 
equal-time moments and the application of the quantum correlation theorem is sufficient 
for the computation of this kind of correlations. 

The expected signal rate in a given experiment averaged over a finite range of energy 
and time can be described in general by a simplified expression of the form: 

N = [ dEdt <j){E) P(E;t) (32) 
JAE,At 

The function <f>(E) is the experimental effective flux, where the original neutrino flux, 
the experimental detection and geometrical efficiencies and cross sections are included. A 
hidden dependence, potentially relevant, on Eo,to, the central points of the integration 
intervals, is understood. 

For a fluctuating survival probability P(E;t), N is itself a random ,non Gaussian, 
variable. The averaged signal is simply the integral of the averaged probability. The 
second moment of the N distribution or more complicated energy or time correlations, i.e. 
averages of the type 

(N(E' ,t' )N(E ,to)) 
are all of the same form, schematically: 

(N 2 ) = [ dE h2 dt l ^{E l )ct>(E 2 ){P(E l ,t l )P{E 2 ,t 2 )) (33) 
J AE,At 



The correlation function appearing in the integrand of Eq.(|33|) can be computed in 
two steps. First we define the generalized density matrix as the tensorial product of 
usual density matrices at two different times and energies: 

X^(E 1 ,E 2 ;t 1 ,t 2 ) = X(E 1 ,t 1 )^X(E 2 ,t 2 ). (34) 

(2) 

The average of the element (Aj^) = (An An) is evidently the probability correlation 
function we are looking for. 

The differential equation for the equal time X^ (Ei,E 2 ;t, t) is obtained from the indi- 
vidual evolution equations for the matrices X\ >2 = X{E\p) ( indices (1,2) label respectively 
expressions where Ei,E 2 appear) 

d t X^ = HiXW +xWH 2 = pLXW (35) 

Eq.(^) is a random linear differential equation, linear in the stochastic variable p(t). 
Applying the formalism developed in the previous section, we can write immediately the 
equation for the ensemble average (A^ 2 )} (Eq.(|lS-14)). Once we know the equal time 



correlator {P(t)P(t)) 7 we obtain the expression for any other pair (t, t') using the quantum 
regression theorem (|l6|]) which reads as follows. For any vector Markov process Y, if the 
ensemble average of Y fulfills a Schroedinger-like equation of the type 

d t (Y(t)) = G(t)(Y(t)), (36) 
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with G(t) an arbitrary matrix, then the second order correlations will obey the following 
equation 

d T (Y i (t + T)Y l (t))=J2G ij (r)(Y j (t + T)Y l (t)). (37) 
j 

Note that the quantum regression theorem is in principle not applicable to systems 
described by the integral equation Eq. (|i~3l) . Nevertheless it is applicable to problems where 
the integral equation can be reduced to a ordinary differential equation as Eq. (|3l|) . Such 
second order equation can be easily expressed as a first order one by defining the auxiliary 
pair process Y = (X, dtX). 

For the simplest case where p is a 5 correlated process the equation for non-equal time 
correlators is explicitly given by: 

d T {X {2 \t + T,t)) = -iQ 2 L 2 (T)(X^(t + T,t)}. (38) 



For exponentially correlated problems second order equations similar to Eq.(31) can be 
immediately obtained. 



7 Conclusions and final Remarks. 

In conclusion, in this work we have developed the general formalism for the study of 
neutrino propagation in stochastic media. This formalism has allowed us to derive sys- 
tematically evolution equations for averaged quantities as survival probabilities and higher 
order moments of them. New integro-differential equations valid for finite correlation pro- 
cesses have been obtained for the first time. For exponentially correlated processes a 
second order ordinary equation is obtained as a consequence. The Redfield equation valid 
for Gaussian delta-correlated noise is rederived in a simple way as a particular case. In 
the context for neutrino oscillations, it has been shown that the presence of matter noise 
induces the appearance of a effective complex part in the matter density. The presence 
of magnetic field noise affects in a qualitatively different way. The wide range of validity 
of the 5-correlated model has been checked in a realistic example corresponding neutrino 
oscillations in random solar matter. Finally, the formalism together with the quantum 
correlation theorem has been used to obtain equations for higher distribution moments 
and non-equal time correlation functions. 

The formalism can be generalized in an obvious way to perform the ensemble average 
of equations of slightly more general type than Eq.(|l|), equations of the form: 

id t X = ^ A (t)L,(t)X; X(0) = X . (39) 

i 

These equations appear for example in chaotic neutrino magnetic precession. The result 
for the case where the pi are (^-correlated in time but mutually uncorrelated is simply: 

id t (X(t)) = -iJ2^L 2 (t)(X(t)). (40) 

i 
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The generalization to continuous Schroedinger equations in presence of random poten- 
tials or random external forces is also obvious if we left aside the mathematical differences 
coming from the appearance of infinite dimensional Hilbert spaces. The formalism is of 
general application to any quantum system governed by similar Schroedinger equations: 
equations where a stochastic function appears multiplicatively in some term of the Hamil- 
tonian. 
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(A) 




Figure 1: Neutrino oscillations in presence of random solar matter. (A) Averaged survival 
probability as a function of Am?/2E ( <5-correlation case, Eq.(|n| )). Continuos curve: 
absence of noise. Dashed lines; f2 2 = 100, 1,0.1,0.01 (x 1000 Km, from longer to shorter 
dashing). (B) The neutrino survival probability is plotted as a function of L/d with 
L = 1/e, "d" is the distance over which the chaotic field is acting (~ 0.6 Rq). Eq.(|3l"|) is 
used for the same values of f2 2 as before ( Am 2 /2E = 6.3 10 -8 eV 2 /MeV now). For both 
plots, the neutrino, created at r « 0.4 R Q as a pure electron neutrino, propagates radially 
outwards. The flavour mixing angle is sin 2 26 = 0.05. 
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